Percolating through networks of random thresholds: 
Finite temperature electron tunneling in metal nanocrystal arrays 
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We investigate how temperature affects transport through large networks of nonlinear conductances 
with distributed thresholds. In monolayers of weakly-coupled gold nanocrystals, quenched charge 
disorder produces a range of local thresholds for the onset of electron tunneling. Our measure- 
ments delineate two regimes separated by a cross-over temperature T* . Up to T* the nonlinear 
zero-temperature shape of the current-voltage curves survives, but with a threshold voltage for 
conduction that decreases linearly with temperature. Above T* the threshold vanishes and the low- 
bias conductance increases rapidly with temperature. We develop a model that accounts for these 
findings and predicts T* . 
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Densely packed arrays of nanocrystals form a new class 
of "artificial solids" with tunable electronic, magnetic 
and optical properties [1-7] . These properties stem from 
single-electron charging and quantum confinement ener- 
gies on the individual particle level, mediated by the cou- 
pling to neighboring particles [8]. However, even for the 
simplest such artificial solid, a close-packed monolayer of 
metallic nanoparticles, the effect of temperature on the 
overall transport properties is poorly understood. 

At low temperatures, metal nanoparticle arrays show 
pronounced nonlinear current-voltage (IV) characteris- 
tics, with a well-defined threshold above which the ap- 
plied bias has to be raised for conduction to occur. Sim- 
ilar threshold behavior occurs in a large family of disor- 
dered nonlinear systems, including the depinning of mag- 
netic flux bundles in type-II superconductors [9] and the 
onset of motion in charge-density waves trapped by im- 
purities [10]. In metal nanoparticle arrays the disorder 
arises, even for perfectly ordered particle arrangements, 
from quenched background charges in the substrate or 
the matrix surrounding the particles. These quenched 
charges induce random shifts in the local nanoparticle po- 
tentials and lead to a distribution of electrostatic energy 
costs incurred by electrons as they tunnel from particle to 
particle. The system behaves as a network of ultrasmall, 
Coulomb-blockade-type tunnel junctions [11] with a dis- 
tribution of threshold energies. For T = 0, a key result 
was obtained a decade ago by Middleton and Wingreen 
(MW) [12] who showed through simulations and a scaling 
model that the effective voltage threshold V t is propor- 
tional to the array length and that the current-voltage 
(IV) characteristics exhibit power law scaling of the total 
current, I, with excess voltage, V — Vt- These predictions 
are in line with experimental results on a wide range of 
nanocrystal arrays [1,4,5,7,8]. More recently, the effect 
of structural disorder on the scaling behavior was inves- 
tigated by experiments and reproduced by simulations 



[7,13,14]. 

For finite temperatures, on the other hand, the lit- 
erature to date shows no consensus. In single-particle 
systems the shape of the IV curves near threshold is 
exponentially sensitive to temperature [11]. For ar- 
rays, Heath's group reported that increasing temperature 
rapidly rounds out any nonlinearities [5,8,15]. Black et 
al. also observed linear behavior above 40K [1]. Other 
recent work found that finite temperature had a much 
weaker effect on the shape of the IV characteristics and 
on the threshold V t [7,14,16,17]. Theoretical approaches 
that could calculate temperature-dependent IVs for large 
arrays are currently not available. 

To address this issue, we have performed systematic 
experiments on large, well-characterized two-dimensional 
arrays of gold nanocrystals. Our results show that the IV 
curves retain their low-T, nonlinear scaling properties to 
remarkably high temperatures and that Vt decreases only 
linearly with T. We can explain these findings by assum- 
ing that there is a subset of junctions, percolating across 
the array, for which thermal fluctuations effectively have 
removed the Coulomb blockade. The starting point is 
the picture developed by MW of charge fronts propa- 
gating through the array [12]. However, our model goes 
beyond MW in three key aspects: it introduces a way 
to deal with the effect of temperature on a distribution 
of charging energies, it goes beyond square lattices and 
applies to arbitrary regular networks, and it can explic- 
itly take into account capacitive coupling between neigh- 
boring islands, as is appropriate for close-packed arrays. 
We believe that this very general approach may also be 
applicable to other systems with distributions of local 
threshold values. 

Our arrays consisted of monolayers of 1-dodecanethiol- 
ligated gold nanocrystals, synthesized and deposited 
from solution onto silicon nitride "window" substrates 
[18,19]. Several arrays were deposited onto substrates 
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with thicker SiaN4 layers and no windows, more suit- 
able for measurements above 200K. Eleven arrays were 
studied ranging from monolayers with long-range hexag- 
onal order, with domains of typical size equal to the elec- 
trode spacing, to monolayers with some disorder, includ- 
ing voids or localized regions of double layers not more 
than 20% in area fraction. In each sample, the nanocrys- 
tals were monodisperse to within 5% as checked by trans- 
mission electron microscopy (TEM). Different samples 
had mean core diameters 4.5nm < d < 6.5nm, and mean 
spacings 1.5nm < s < 2.6nm. The relevant Coulomb 
charging energy, AE, is the change in electrostatic en- 
ergy of the system due to an electron tunneling from one 
particle to the next. AE depends on the background 
charges on all particles involved. However, because of 
compensation by mobile electrons, it suffices to consider 
effective background charges within [-e/2, +e/2] and the 
distribution of charging energies, P(AE), only extends to 
±AE max . For close-packed monolayers, keeping only in- 
teractions between the two particles forming the junction 
and the surrounding 8 nearest neighbors, to good approx- 
imation AE max = e 2 (C^ + 2C 12 1 ) = e 2 Cf 1 1 (l + 2 7 ), 
with 7 = Cf 2 /Cn where C _1 is the inverse capaci- 
tance matrix. From measurements of d and s, and tak- 
ing a dielectric constant e = 2 for the alkanethiol, the 
inverse capacitance matrix elements for such 10-particle 
system were calculated numerically [20]. This gave Cj^ 1 
of order 10 18 F _1 and typical values for 7 « 0.36. All 
arrays had a fixed width (2/im, or w 270 particles) and 
lengths, N, ranging from 27 - 170 particles, set by the 
electrode separation. For samples on window substrates, 
N was measured by direct inspection with TEM after 
the transport measurements were completed. Otherwise, 
N was estimated from the electrode separation. Samples 
were cooled in vacuum, and dc IV curves were taken 
using Keithley 614 electrometers at voltage ramp rates 
5 - 25mV/s. 

At low temperatures, all arrays showed a clear volt- 
age threshold for conduction, indicating strong Coulomb 
blockade behavior. Close inspection of the TV curves 
(see Fig.l) reveals that the differential conductance, g = 
dl/dV, for a given current I is independent of temper- 
ature below a sample-dependent characteristic tempera- 
ture in the range 100 K - 200 K. The TV curves keep 
their shape and simply shift to lower voltages with in- 
creasing temperature. More quantitatively, we find that 
the IVs can be collapsed on top of each other by trans- 
lating them by an amount V s hat(T) < 0. Defining an 
effective threshold V t (T) = vf(0) + V shift {T), we can col- 
lapse all traces onto a master curve as shown in Fig. lb. 
|14hift(T'*)| — V t (0) defines the temperature, T* , at which 
the array threshold first reaches zero. Once T > T* , the 
IVs exhibit non-zero differential conductance at zero bias 
(Fig. la inset) prohibiting a data collapse at small bias. 
However, at higher bias the IVs still collapse when trans- 
lated by an amount |V r s hift(? 1 )| > Vt(0). The regime above 
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FIG. 1. a) Evolution of IV curves with temperature for an 
array of length N — 128. The IVs are fully symmetric with 
respect to bias voltage reversal. Above T* the IVs no longer 
have a finite voltage threshold, as shown in the inset, b) 
Below T* , IVs at different temperatures can be collapsed by 
translating the voltage scale by an amount V s hHt(T), with no 
alteration of the current scale. The inset shows the zero-bias 
conductance, go = dl/dV, as a function of temperature. 

T* thus is characterized by a negative effective threshold 
and low-bias data falling below the master curve. 

For highly ordered arrays (bottom curve in Fig. 2a) 
MW's T = scaling law, I ~ (V - V t ) c , remains 
valid up to T* with temperature-independent exponent 
C = 2.25 ± 0.1. Details of the functional form of I(V) 
in different arrays affect the value of T* but are irrele- 
vant as far as the collapse is concerned. Even for slightly 
disordered arrays, characterized by a somewhat s-shaped 
IV [7] in Fig. 2, collapse can be achieved. For all arrays 
studied, V t (T), the only adjustable parameter in the scal- 
ing procedure, drops essentially linearly with increasing 
temperature (Fig. 2b). Both 14(0) and the slope of this 
depression increase with array size N. Furthermore, for 
similar V t (0) the slope decreases with increasing charging 
energy. 

This slow, linear depression of the voltage threshold 
with temperature is surprising and unexpected from the 
theory of few-junction networks [11] or from other nonlin- 
ear systems such as charge-density waves [10]. A similar 
trend was first noted in ID chains of carbon nanospheres 
by Bezryadin et al. [16], and by Ancona et al. [14] in 
small 2D arrays. Bezryadin et al. proposed that ther- 
mal fluctuations reduce the charging energy by kgT in 
all particles and that therefore V t (T) should decrease as 
NkgT/e. However, this simple picture cannot account 
for the slope values in Fig. 2b. We now outline a model 
that captures the linear dependence of V t (T), predicts 
the slope, and, under appropriate rescaling, collapses the 
data in Fig. 2 onto a universal curve. 

MW showed that the minimum bias voltage for charges 
to move all the way from one electrode to the other de- 
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FIG. 2. a) Scaling behavior of current I with excess voltage 
V — V t (T) above threshold for same array as in Fig.l (middle 
curve) and two other arrays with different degrees of disorder. 
For each array, 4 different temperatures are plotted. The 
dotted line corresponds to ( = 2.27. b) The effective threshold 
as a function of temperature for several arrays of length N and 
effective charging energy eVt(0)/aN. 

pends on the number of "up-steps" along any given path. 
Each local up-step in potential, no matter how small, re- 
quires an increment AV in the externally applied bias for 
the charge front to advance, and thus contributes equally. 
(AV = AE max /e for uncoupled particles or junctions 
surrounded by charge neutral neighbors; for the more 
general case see below.) The array threshold is set by 
the optimal path with the smallest total number of up- 
steps and scales as Vt — aN AV in both ID and 2D [12]. 
In 2D, a = 0.338 for square lattices [12] and a = 0.227 for 
triangular lattices [21]. Down-steps do not gain energy 
because all tunneling between particles is inelastic (no 
resonant tunneling between sharp levels). Hence, as long 
as the number of up-steps along a path remains fixed, 
a gradual reduction in the associated energy step size, 
induced by thermal fluctuations, will not reduce Vt- We 
now argue that V t can be reduced if temperature elim- 
inates some of the up-steps completely. This occurs for 
the subset of junctions with step magnitudes \AE\ so 
small that they are washed out by thermal fluctuations. 

We consider a scenario in which temperature has re- 
moved the nonlinear, threshold-type Coulomb blockade 
behavior of a fraction p(T) of the junctions (irrespective 
of the sign of the step), effectively turning them into lin- 
ear resistors. The effect of temperature, then, is a shift 
in threshold proportional to p(T). At some temperature, 
T*, the fraction p(T*) will become large enough for a con- 
tinuous path of junctions without discernable Coulomb 
blockade to form a percolating sub-network that bridges 
the array. At this point the overall array voltage thresh- 
old will vanish. We therefore identify p(T*) with the per- 
colation threshold, p c , for the underlying lattice. Hence 
we expect V t (T) = V t (0)(l - p(T)/p c ). For a triangular 



lattice p c « 0.347. To make explicit the temperature de- 
pendence of the threshold voltage, we need to find p(T). 
For a given distribution, P(AE), junctions corresponding 
to \AE\ < bk-sT are taken as effectively linearized. The 
factor b accounts for the broadening of the Fermi-Dirac 
distributions in the nanocrystals. Electrons in the high- 
energy tail of one nanocrystal's distribution can tunnel 
into the low-energy tail of accessible states of a neighbor- 
ing particle, reducing the energy cost by roughly 2k^T. 
More accurately, equating the reduction with the differ- 
ence between the mean energies in the high- and low- 
energy tails, we calculate b = 2.4 [21]. At any given tem- 
perature, p(T) is then given by the integral over P(AE) 
between the limits ±bksT. 

Consider first the case of two particles with random 
charge offsets in the range [— e/2, e/2], forming a repre- 
sentative tunnel junction, surrounded by 8 neutral neigh- 
bors. In this case, AE max = e 2 C^(l - 7), and P(AE) 
is a triangle with y-intercept P(0) = l/AE max (inset 
to Fig. 3) [21]. Integrating P(AE) over the shaded area 
leads to p(T) = bx(2 - bx) with x = fc B TP(0). Us- 
ing p(T*) = p c and b = 2.4 we find k B T* « 0.08/P(0). 
Note that p(T) is proportional to T and -P(O) as long 
as x << 1, which is satisfied up to temperatures far 
above T*. Allowing for random offset charges also on 
the 8 neighbors rounds the corners and steepens the 
slope of P(AE). Still, the overall triangular shape re- 
mains, p(T) = 2bx + C(x 3 ) and in the above consid- 
erations only the numerical value for P(0) needs to be 
modified. Now AE max = e 2 C 1 " 1 1 (l + 2 7 ) and P(0) = 
(e 2 Cf 1 1 )- 1 (l - 1.57 7 )/(1 - 7 ) 2 [21]. For all arrays, inde- 
pendent of N, the model then predicts that the normal- 
ized threshold Vt(T)/Vt(0) decreases linearly as a func- 
tion of the dimensionless variable x — ksTP(0), with a 
slope —2b/p Cl crossing zero near x = 0.072. Within this 
picture, 1 /P(0) can be thought of as an effective charging 
energy per up-step that captures all effects of the local 
packing geometry and the quenched charge distribution. 

These predictions agree well with our data (Fig. 3) [22]. 
To obtain V t (0) we extrapolated V t (T) back to T = 
for each sample (Fig. 2b). To determine P(0) two meth- 
ods were used, giving independent checks of the model. 
Equating eAV = 1/P(0) for the energy per up-step we 
obtained P(0) from Vt(0) = aN AV using N from the 
electrode separations and a = 0.227. In addition, for 
the samples where TEM images of sufficient resolution 
were available to extract d and s, we were able to cross- 
check the values for P(0) with those derived from direct 
numerical calculation of the C" 1 [20]. 

The broad distribution of local energy differences AE 
guarantees that, even if a fraction of steps is effectively 
removed, the underlying concept of current paths me- 
andering through a random potential landscape - the 
premise of MW's scaling picture - remains valid. The 
general shape of the IVs above V t (T), which reflects how 
additional current paths open up with increasing bias be- 
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FIG. 3. Normalized threshold, V t (T)/V t (Q) as a function 
of dimensionless temperature variable x = fcsTP(O) for all 
arrays. The dashed line is the model prediction. Inset: Prob- 
ability density of charging energies P(AE). With capacitive 
coupling in the presence of offset charges on all nearest neigh- 
bors, P(AE) evolves to a more rounded shape. 

yond threshold, therefore should remain unchanged up 
to T* . This is borne out by our data. Rounding of the 
Coulomb blockade in the local junction IV characteris- 
tics affects the quality of the data collapse only at the 
smallest V - V t (T) (Fig. 2a). 

The behavior for T > T* is quite different. Here the 
low-bias conductance is dominated by a few junctions 
from the (p c /2)-th percentile of P(AE), since they are 
the key, "bottleneck" bonds in the lattice during perco- 
lation at T*. Therefore, the low-bias IVs for the perco- 
lating path resemble those of a single island system with 
washed-out Coulomb blockade, for which the linear zero- 
bias conductance (Fig. la inset) exhibits thermally acti- 
vated behavior [11]: go(T) ~ cxp (-U/k B T). Our model 
predicts U = bk B T* = Q.17AE max , the largest AE for 
the percolating path. This g will dominate the low-bias 
array conductance as T increases and additional, parallel 
paths open up that have larger U. Because of these addi- 
tional paths, we expect the measured effective activation 
energy to be somewhat larger than bksT* . Our data in- 
deed are compatible with simple activated behavior and 
a magnitude of U 1 to 2 times bksT* . For example, 
for the array in Fig. lb inset AE max = 124 meV, giving 
bksT* = 21 meV, to be compared with U w 38 mcV 
from the slope of In (go) vs T _1 . However, because of 
the limited (high) temperature range available, our data 
cannot rule out other functional forms, such as variable 
range hopping conduction as seen in Ref. [15]. 

Our results support a picture of thermally assisted 
quantum tunneling between neighboring gold nanocrys- 
tal cores through the dodecanethiol ligand. Within this 
picture, the ligands serve as a mechanical spacer between 
the particles and they set the tunneling barrier height. 



They do not introduce additional states into the barrier, 
nor is the barrier low enough to be thermally hopped 
over [23] , both of which may be the case with more com- 
plex ligands [4]. We find no evidence of a regime with 
positive temperature coefficient of resistance as seen in 
Rcf. [5]. Such behavior may require a stronger inter- 
particle coupling and junction resistances much closer to 
the quantum resistance R q = h/e 2 . In the weak coupling 
limit, as in our arrays, the value of T* as well as the col- 
lapse of the temperature-dependent threshold data onto 
a universal curve emerge naturally from considering the 
interplay of tunneling and quenched charge disorder. 
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